Diffusion of two molecular species in a crowded environment: theory and experiments 
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Diffusion of a two component fluid is studied in the framework of differential equations, but where 
these equations are systematically derived from a well-defined microscopic model. The model has a 
finite carrying capacity imposed upon it at the mesoscopic level and this is shown to lead to non-linear 
cross diffusion terms that modify the conventional Fickean picture. After reviewing the derivation 
of the model, the experiments carried out to test the model are described. It is found that it can 
adequately explain the dynamics of two dense ink drops simultaneously evolving in a container filled 
with water. The experiment shows that molecular crowding results in the formation of a dynamical 
barrier that prevents the mixing of the drops. This phenomenon is successfully captured by the 
model. This suggests that the proposed model can be justifiably viewed as a generalization of 
standard diffusion to a multispecies setting, where crowding and steric interferences are taken into 
account . 

PACS numbers: 05.60.Cd, 05.40.-a, 87.15.Vv 



I. INTRODUCTION 

Diffusion is a key process in nature which describes 
the spread of particles from regions of high density to 
regions of low density. As is usually the case, the first 
quantitative approach to the study of mass transport via 
diffusion was phenomenological, by Fick in the middle of 
the 19 th century. His first law relates the diffusive flux 
J(x, t) to the concentration 0(x, t) at position x and time 
t, and mathematically amounts to the relation 



-DVt 



(1) 



where D is a constant. Equation ([T}, together with mass 
conservation, yields a diffusion equation for the concen- 
tration 6: 



1t= DV 



(2) 



which is also referred to as Fick's second law 

A microscopic justification for the diffusion equa- 
tion followed some time later: Einstein Q and Smolu- 
chowski Q independently proposed a physical interpre- 
tation of the phenomenon, although based on the experi- 
mental work of Brown Q almost a century earlier. They 
showed how the diffusion constant, D, can be related to 
measurable physical quantities such as viscosity and tem- 
perature. In addition, the famous relation between the 
mean square displacement ((x 2 )} and the time elapsed t, 
((x )) = 2Dt, can be derived. Here double angle brack- 
ets indicate cumulants. The unintuitive, but now univer- 
sally accepted, conclusion is that the average root-mean- 
square displacement of a Brownian particle grows with 
the square root of time. 



Despite the fact that the above picture, and Fickean 
diffusion in particular, is generally appropriate when de- 
scribing the spontaneous spatial rearrangement of parti- 
cles in suspension, deviations are expected to occur in 
various situations of interest e -S- if obstacles are 

present or when different, and so distinguishable, species 
are sharing the available space at high concentration. A 
vast literature reports on these anomalies (see Q and ref- 
erences therein); they are particularly important in ap- 
plications to molecular biology, where a large ensemble of 
microscopic entities (e.g. proteins, macromolecules) are 
densely packed in space. This crowding effect can for ex- 
ample force the molecules populating the cells to behave 
in radically different ways than in test-tube assays 0- 
[Tl| . For this reason, the study of biochemical processes 
under realistically crowded conditions is absolutely cen- 
tral to the understanding of cellular environments. These 
operating conditions are in fact ubiquitous in cells, and 
crowding could be an essential ingredient for an efficient 
implementation of the metabolic machinery. 

When crowding occurs — a macroscopic effect reflect- 
ing the microscopic competition for the available spatial 
resources — the mean-square displacement of particles is 
seen to scale in a sub-diffusive fashion. This empirical 
fact has so far being described by resorting to different 
fitting strategies for the experimental data. The most 
popular of these is to hypothesize a power law scaling for 
the variances [l2[. Fractional diffusion equations have 
also been proposed to explain the emergence of the scal- 
ing invoked [7|, but the approach frequently appears to 
be a mathematically sound expedient which reproduces 
the observed behavior, but without a firm physical basis. 

Interestingly, practically all of the life science appli- 
cations focusing on the formation of spatial and tempo- 
ral patterns, assume reaction-diffusion schemes inspired 



2 



by the mathematical paradigms of single species diffu- 
sion. As such, the underlying partial differential equa- 
tions have only a simple Laplacian term, and no accom- 
modation is made for nonlinear corrections and/or cou- 
plings that need to be introduced if obstacles, or more 
than one species, are present at hi gh c oncentration. 

By contrast, in an earlier work [131 ] a modification of 
Fick's law was proposed which was based on a detailed 
microscopic theory, which effectively encapsulated the es- 
sential aspects of limited spatial resources. This first- 
principles procedure represents a systematic approach to 
the modeling of biological, biochemical and related sys- 
tems; it begins from an individual-based description of 
the interaction between the microscopic elements which 
make up the system. The competition for the available 
resources leads to a modified (deterministic) diffusive be- 
havior. More specifically, cross-diffusive terms appear 
which link multiple diffusing components and which mod- 
ify the standard Laplacian term, a relic of Fick's law. The 
finite carrying capacity of each microscopic spatial patch 
is the physically motivated assumption at the heart of 
this extended theory of diffusion. The theory has been 
used in the study of reaction-diffusion systems which are 
relevant to applications in the life sciences [3, [H| • Other 
models have been proposed which, starting from micro- 
scopic rules, lead to modified macroscopic diffusion equa- 
tions fill ]. 

This paper discusses an experimental investigation of 
the generalized model of diffusion presented in [l]| . To 
this end, we will study the diffusion of ink drops lo- 
cated in the same spatial reservoir. As we will show, 
the diffusion of an isolated drop is successfully described 
in terms of standard diffusion. The case of different types 
of ink drops simultaneously evolving is more difficult to 
describe: the diffusion comes to an halt when contiguous 
drops make contact. This regime, clearly non-diffusive, 
is explained by invoking the multispecies setting of (l3j— 
an observation that we shall quantitatively substantiate 
by benchmarking measurements to simulations. 

The paper is organized as follows. The next section, 
Section II, is devoted to introducing the model and re- 
viewing the main steps of the derivation given in [l3|. 
In Section III we will discuss the experimental apparatus 
and, in Section IV focus on the evolution of a single ink 
drop inside a Petri dish containing water. The coevolu- 
tion of two drops is studied in Section V, where fits to 
the data allow the accuracy of the theory to be assessed. 
Finally, in Section VI we will sum up and conclude. 



II. THE MODEL 

In the following we shall review the derivation of 
the model in the simplified case where two species of 
molecules diffuse, sharing the same spatial reservoir. It 
is however worth emphasising that the analysis extends 
straightforwardly to the generalized setting where an ar- 
bitrary number of molecular species are present. Sim- 



ilarly, the analysis holds in an arbitrary number of di- 
mensions, but below we will assume that we are working 
in the physically relevant case of two dimensions. 

The generic microscopic system that we consider will 
then be constrained to occupy a given area of two- 
dimensional space. We assume that this area is parti- 
tioned into a large number f2 of small patches or cells. 
Each of these mesoscopic cells, labeled by i, is charac- 
terized by a finite carrying capacity: it can host up to 
N particles, namely n\ of type A^ l \ of type A^ 2 \ 

and Vi = N — n- 1 ' — nf^ vacancies, hereafter denoted by 
V. The molecules are assumed to have no direct inter- 
action. There is however an indirect interaction which 
arises from the competition for the available spatial re- 
sources. The mobility of the molecules will be de facto 
impeded if the neighboring cells have no vacancies. More 
concretely, we shall assume that molecules move only to 
nearest-neighbor cells, and only if there is a vacancy that 
can be eventually filled. This mechanism translates into 
the following chemical equation 
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(a) 



.(a) 



V t + A 



1,2, 



(3) 



where i and j label nearest-neighbor cells with , A\ 2 ^ , 
and Vi identifying the particles (including the vacancies) 
in cell i. The parameters /u,^ and fi^ denote the as- 
sociated reaction rates. The state of the system is then 
characterized by the number of A^ and A^ molecules 
in each cell; the number of vacancies obeying a trivial 
normalization condition. To make the notation com- 
pact, we introduce the vector n = (ni, . . . , tiq), where 

rij = (n^,n^'). The quantity T(n'|n) represents the 
transition rate from state n, to another state n', com- 
patible with the former. The transition rates associated 
with the migration between nearest-neighbors, as con- 
trolled by the chemical reaction Eq. ([3]), can be cast in 
the form 



T(n^- l,nf> + l\nf>,nf>) 



» 



( a ) „( a )) 



zQ, N N 



= 1,2, 
(4) 

where z is the number of nearest-neighbors that each cell 
has. In addition, and to shorten the formulae, in the 
transition rates we only show the dependence on those 
particles which are involved in the reaction. As discussed 



in [13[ , it is the factor Vj = N — n 



(i) 



(2) 



reflecting the 



finite carrying capacity of cell j, that will eventually lead 
to the macroscopic modification of Fick's law of diffusion. 

We assume that the process is Markovian, so that the 
probability of the system being in state n at time t, 
P(n, t), is given by the master equation [l7| 

^JM= £ [T(n|n')P(n' ) t)-T(n'|n)P(n ) t)] ) 

n'-^n 

(5) 

where the allowed transitions are specified by Eq. (QJ. 
Starting from this microscopic, and hence inherently 
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stochastic picture, one can recover the corresponding 
macroscopic, deterministic description. To do this we 
need to determine the dynamical equations for the en- 



semble averages and (n^). This is achieved by 

multiplying the master equation by nj a ^ and sum- 
ming over all n. After an algebraic manipulation which 
requires shifting some of the sums by ±1, one finds 



dt 



Y:[(T(n^ + l,n^-l\n^,n^)) 



jet 



(T(n[ a) -l,nf + l\n\ a ',nf')) 



(°) „,( a h 



(6) 



where the notation Yljei means that we are summing 
over all patches j which are nearest- neighbors of patch i. 

The averages in Eq. (JSJ are explicitly carried out by 
recalling the expression for the transition rates (|4|) and 
then replacing the averages of products by the products 
of averages, a formal step which is legitimate in the con- 
tinuum limit N — > oo [17| . Rescaling time by a factor 
Nfl, gives [Hi 



dt 
dt 



,(!) 



,(2) 



A<A 



A0 



(7) 



In the above expression we have introduced the con- 



centration 



= limjv->oo js ■ The symbol A de- 
notes the discrete Laplacian operator which reads Afi = 
(2/z) Yljeiifj ~ fi)- Finally, by taking the size of the 
cells to zero, and scaling the rates // a ' appropriately [HI 
to obtain the diffusion constants D' a ', one finds the fol- 
lowing partial differential equations for the continuum 
concentration ^^(x, t) and 4>( 2 \x,t): 



dt 



dt 



72 J, 2) 



(2) vV 1 ) -4,W v 2 d>W 



(8) 



where V 2 is the usual Laplacian. 

The extra terms ±(0^) V 2 ^ -^V 2 ^ 1 )) in the dif- 
fusion equations ([5]) arise directly from having included 
vacancies in the transition probabilities. These nonlin- 
ear modifications of the simple diffusion equation have 
therefore a specific microscopic origin, as opposed to 
various ad hoc suggestions that have sometimes been 
proposed. These equations cannot be solved exactly; 
we have studied them numerically using both the for- 
ward Euler method in time and a semi-implicit algorithm 
which exactly conserves the number of particles [19] , and 
checked that they agree. The explicit results given in this 
paper were obtained using the forward Euler method. 

Since the two types of molecules are separately con- 
served, and </>( 2 ) satisfy the conservation relations 



dcj^Jdt + divJ^(i) = and d^jdt + divJ 0(2) = 0. 
This means that the fluxes Jmi) and Jm2) satisfy 



The above relation generalizes Fick's first law (p}. 

In [l3j a quite complete analysis of Eqs. © is given in 
the simplified case when the two diffusion constants are 
equal. It is useful to define $( a) = / (a) dx, a = 1,2, 
where quantifies the amount of material of type 

. The integrations are carried over the space where 
the co-evolving molecular species are confined. In the 
simulations performed later, this will be a square do- 
main of size L. It is therefore useful to introduce 0( a ' = 
<$>^/L 2 . Here, 4 a > can be interpreted as a type of (non- 
dimensional) average number density. 

In this paper, our aim is to provide an experimental 
verification of the generalized diffusion equations as ob- 
tained above. To reach this goal, we have constructed a 
simple experimental apparatus that enables us to quan- 
titatively track the evolution of ink drops in water, and 
specifically considered the case of interest where two 
drops arc simultaneously present. The theory will be 
tested against direct measurements so as to assess its va- 
lidity. The next section is devoted to introducing the 
experimental set-up which we use. 



III. THE EXPERIMENTAL APPARATUS 

The experiments that we will discuss in this paper are 
aimed at monitoring the progressive diffusion of ink drops 
when deposited in water. This has be en p reviously stud- 
ied systemically only for a single drop [20j . We will repro- 
duce these results, but only as a prelude to the real case 
of interest: the interaction between two drops. In this 
Section we will briefly describe the experimental appara- 
tus, as well as the procedure of data acquisition which 
we used. 

A Petri dish of diameter 14cm was filled with a uniform 
layer (height 0.5cm) of water. The water was treated with 
a small amount of surfactant, to reduce the effects of sur- 
face tension. The container was placed on a horizontal 
surface, illuminated from below with a light emitting ca- 
pacitor (LEC) produced by Ceelite [2l[ and distributed 
by Continua Light [22J . This was the only source of light 
employed during the course of the experiment. Using 
an illumination in transmission strategy enabled us to 
substantially reduce the impact of spurious, undesirable 
shadows on the recorded images. The LEC panels were 
lightweight and extremely flexible; those in the experi- 
ment had dimensions of 21cmx35cm. A color CCD cam- 
era (DFK 21AF04) [H] was placed 40cm above the con- 
tainer and controlled remotely via the computer, with an 
image of the sample being taken every 5 seconds. A pho- 
tograph of the experimental set-up is displayed in Fig. [T] 
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FIG. 1: (Color online) The experimental set up. The cam- 
era is suspended above the Petri dish and the system is il- 
luminated from below by a light emitting capacitor. The si- 
multaneous insertion of the two drops is carried out with a 
mechanical rod which holds two graduated injectors together. 
This device allows for the change of the distance between the 
drops at will, and in a completely reproducible fashion. 

The solution was left to settle for several minutes so as 
to damp out any source of uncontrolled perturbations, 
such as local vortices or turbulent flows, that may de- 
velop when the water is put into the dish. Single, or 
multiple, drops of ink were then introduced with a grad- 
uated injector. The ink used in the experiments was of 
the Ecoline type, Talens brand [24 |. Two specific col- 
ors were used: magenta and cyan. Upon injection, the 
ink formed a roundish drop which, under the action of 
gravity, fell to the bottom of the Petri dish. The system 
was therefore essentially two dimensional, the thickness 
of the ink layer being very small (of the order of tens 
of microns). The drop now started growing: the circu- 
lar profile being preserved for isolated drops, and, in this 
case, the average drop radius progressively increased with 
time. By contrast, when two or more drops were present, 
mutual interference occurred that slowed down the diffu- 
sion, which was translated into distorted profiles. These 
latter profiles are analyzed in Sec. V in the simplified 
case where just two drops were made to simultaneously 
evolve. In the next section we report on the experiments 
carried out on isolated drops. 



IV. DIFFUSION OF ONE ISOLATED DROP 

The simplest situation, and the one which we used to 
calibrate parameters associated with the ink, was when 
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FIG. 2: (Color online) Magenta drop: a gallery of images 
(a)-(t) is displayed at different times covering 495s with ac- 
quisition step of 5s. The images are treated so as to isolate 
the region of interest where the drop is localized. The exper- 
iments are carried out at a temperature of 17 ± VC. 



one isolated drop was allowed to form within the Petri 
dish [2(|. The ink was injected at one specific location; 
varying the amount of material that was introduced did 
not alter the fact that it formed into a regular, almost two 
dimensional, circular domain. The region of interest in 
the recorded image was then cropped, and the operation 
repeated for all successive snapshots of the dynamics, re- 
sulting in the gallery of images displayed in Fig. [2] A 
movie of the dynamics of the drop is attached as supple- 
mentary material. 

The photographs that were obtained were then pro- 
cessed using the freely available Image J software [25j | . 
Depending on the color of the ink, a specific filter was 
chosen, which maximized the corresponding adsorption 
and so enhanced the contrast of the drop. Thus, the 
green channel was employed when analyzing the magenta 
signal, while the red channel was used for the cyan drop. 
This choice simply follows from the subtractive theory of 
the CMYK color model H^. 

Each image, converted into an intensity matrix, was 
normalized to the value of the maximum and then ex- 
pressed in binary form. This was carried out by choosing 
a given (arbitrary — see later discussion) threshold and 
setting pixels with an intensity lower than this threshold 
to 0, and those with an intensity above this threshold 
to 1. Each binary map was filtered to eliminate isolated 
pixels and small spurious features. 

From this set of binary images, the diameter d m of 
the circular profile corresponding to the magenta drop 
was determined. For definiteness, we shall focus on the 
magenta drop to describe the procedure; the subscript 
m denotes magenta. The data are plotted as magenta 
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FIG. 3: (Color online) The diameter d m of the magenta drop 
plotted versus time. The diameter is expressed in pixels (1 
pixel ~ 0.33mm). Snapshots are taken every 5s. Symbols 
refer to the experimental data, the diameter of the drops be- 
ing calculated via an automatic procedure which transforms 
the two-dimensional intensity plot into a binary matrix. The 
solid line represents the solution of Eq. where the diffu- 
sion coefficient D m and the relative threshold 7 m are adjusted 
via a least square fitting procedure that maximizes the agree- 
ment between theory and experiments. In the simulations, 
at t — 0, the drop is assumed to uniformly occupy a circular 
area. The radius of the circle (in pixels) is assigned so as to 
approximately match the size of the drop as seen in the first of 
the available experimental images. The simulations are per- 
formed over a square domain Lx L, with L — 150 pixels. The 
best fit values are D m = 0.7712 pixels^" 1 and -y m = 0.0675. 



colored circles in Fig. [3J This quantifies the growth of 
the drop, which was qualitatively displayed using the raw 
data in Fig. [5] In Fig. [31 the vertical axis is d m , expressed 
in pixels. The conversion from pixels to physical units is 
straightforward, provided the resolution of the camera 
used is known (see caption of Fig. |3J. 

We hypothesise that the observed dynamics is purely 
diffusive [20| . Motivated by this belief, we seek to estab- 
lish a direct connection between the experimental data 
points, as displayed in Fig.[3J and Eq. ©, which governs 
the dynamics of diffusive processes. It is worth emphasis- 
ing that Eq. © governs the time evolution of the density 
of the diffusing species, while the experimental points 
recorded are calculated from intensity signals as regis- 
tered by the camera. Comparison of the profiles from the 
equation and from the data is therefore not straightfor- 
ward, since the transformation from intensities to densi- 
ties is highly non-linear. It is possible, however, to make 
a comparison at the level of the (circular) domain ob- 
tained from data in the manner described above. The 
idea is to numerically integrate Eq. ^) with the diffusion 
coefficient D (here identified as D m ) being adjusted it- 



crativcly so as to maximize the agreement between the 
theoretical and experimental curves. This is not the only 
free parameter: recall that the threshold imposed when 
performing the binary conversion was arbitrarily chosen. 
Once again, the transformation between that threshold 
parameter and the one which delimits the boundary of 
the drop in the numerical integrations of the equation will 
be highly non-linear. However, we do not need to char- 
acterize this mapping, and need only know that there 
is a one-to-one correspondence between them, using the 
threshold parameter in the numerical integrations as the 
other free parameter of the theory. We denote it by 7 TO . 
Therefore the numerical solution of Eq. ([2]) will contain 
two free parameters: D m and j m . 

We have referred to the solution of Eq. @ as obtained 
from a numerical integration. It is, of course, possible to 
solve Eq. analytically, but we choose to find the solu- 
tion numerically because we wish to use the same method 
in the case of two drops, and to test the methodology in 
the present case of a single drop. Since the correspond- 
ing Eqs. © are not exactly solvable, we are forced in the 
case of two drops to use numerical integration methods. 

The procedure is now that defined above: to iteratively 
adjust the two parameters D m and "f m so as to maximize 
the agreement between the theoretical and experimental 
curves. However, the experimental data gives the radius 
d m of the circular domain, and we need to relate this 
to 7 TO . We choose to do this as follows: if the drop ex- 
tends over a two-dimensional circular domain of radius 
d m , then (f>(x,y,t)/m.ax Xi y((j)(x,y,t)) < 7 rn for each point 
x = (x, y) in the domain. 

The experimental data, represented by circles, are plot- 
ted as function of time in Fig. |3J The solid line represents 
the solution of the diffusion equation @, with the diffu- 
sion constant D m and the threshold parameter j rn found 
using the iterative fitting scheme described above. The 
agreement is excellent, the evolution of the drop being, 
as expected, purely diffusive. Similar conclusions can be 
reached for the cyan drop. Experimental and theoreti- 
cal computed diameters d c are compared in Fig. |4j The 
diffusion constant D c and the threshold 7 C for cyan are 
found by analogy with that for magenta, that is, through 
the iterative fitting procedure described above. 

We conclude that standard diffusion theory adequately 
explains the dynamics of isolated drops, as seen in the ex- 
periments that we carried out. A straightforward fitting 
strategy allowed us to calculate the associated diffusion 
constants for the two different colors of ink that we stud- 
ied. These values will be assumed, and used as input, 
in the more complex scenario where two drops simulta- 
neously evolve in the same Petri dish. We now go on to 
discuss the corresponding analysis in this case. 
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FIG. 4: (Color online) The diameter d c of the cyan drop 
plotted against time. Symbols represent experimental data. 
The solid line represents the solution of Eq. ([2}, where the 
diffusion coefficient D c and the threshold parameter 7 C are 
found via a fitting procedure. The initial condition in the 
simulations is set as described in the caption of Fig. [3] The 
best fit values are D c = 0.9459 pixels^' 1 and 7c = 0.0628. 



V. DIFFUSION OF TWO DROPS: MUTUAL 
INTERFERENCE AND CROSS TERMS 

In this Section we discuss the diffusion of two ink drops, 
one magenta and the other cyan, which were simultane- 
ously inserted into the Petri dish described above. The 
drops were initially localized in different parts of the dish, 
and so did not influence each other. Diffusion made them 
gradually expand, until they eventually came into con- 
tact, as shown in Fig. [5] At this point, a sharp interface 
developed between the drops which acted as a barrier pre- 
venting the mixing of the two compounds. From then on, 
the constraints imposed by this barrier distorted the pre- 
viously circular profiles, as shown in Fig. El We interpret 
this phenomenon as stemming from the competition of 
the molecular constituents for the finite available space. 
The saturation of microscopic spatial resources presum- 
ably results in locally congested configurations, that are 
responsible for the macroscopic features seen in experi- 
ments. 

To test this hypothesis we proceeded as follows. We 
defined an X-axis as the direction of the line connecting 
the centers of mass of the two drops at t = 0. A Y- 
axis was then denned as being orthogonal to the X-axis. 
The images recorded were then processed to measure the 
maximal extension of both of the drops, along the two 
independent directions X and Y , at different times of ac- 
quisition of the signal. This therefore gave four quantities 
at each time t, two for the magenta drop and two for the 
cyan drop. The values obtained were found by filtering 



the recorded intensity map with an imposed threshold, 
following the strategy discussed above. This created a 
binary representation of the two drops, the analysis re- 
turning four time series, two for each drop. The data 
are plotted in Fig. [7] the upper panel displays the data 
corresponding to the magenta drop, the lower panel to 
the cyan drop. In both panels, the circles represent the 
measurements carried out along the Y direction, perpen- 
dicular to the line connecting the centers of mass of the 
drops. The squares refer to the extension of the drops in 
the X direction, along the line connecting the centers of 
mass. 

The initial condition used in solving Eqs. (|5J) mimics 
that found in the experiment. The drops, assumed to 
be approximately circular at t = 0, are in fact assigned a 
specific size that is deduced from the experiments — from 
the first of the recorded images. Similarly, the relative 
separation between the drops is calculated in pixels and 
used as an input to the numerical simulations. 

At relatively short times, before the magenta and cyan 
drops come into contact, their diameters in the X and Y 
directions are approximately the same, their shape being 
approximately circular. As already remarked, once the 
fronts collide (t ~ 100s, see Figs. [5] and [5] panel (m)), the 
dynamics slows down at the interface of the two drops. 
The merging of the drops is consequently impeded along 
the direction of collision, a dynamical effect that con- 
tributes to reducing the growth of the drop width along 
the X direction. By contrast, the diameter of both drops 
keeps on increasing along the Y direction as if normal 
diffusion governed the dynamics. The symmetry of the 
system is therefore broken, as clearly shown in Fig. [7J 

To establish a bridge between the results of the exper- 
iments and the multispecies model of diffusion (|8[). we 
used a fitting protocol based on the analysis developed 
in Sec. IIV1 In particular, the coefficients D^> = D c and 
£)(2) = i) m are obtained from the single drop analysis, 
as are the respective thresholds j c and "f m . Two param- 
eters are left undetermined, the rescaled normalization 
constants and </>( 2 ). These will serve as fitting pa- 
rameters, found by optimizing the agreement between 
theory and measurements. The difference between the 
two, that is minimized by the iterative fitting scheme, 
is the sum of the norm of relative deviations between 
the numerically computed and the experimentally mea- 
sured diameters, along the two directions of observation, 
for both drops, at the same time. It should be empha- 
sised that, due to the nonlinearity introduced by the cross 
terms in Eq. (|SJ), the normalization factors (f)^ and 4>( 2 '> 
do not trivially cancel, as is the case if the linear diffusion 
equation ^ is assumed to govern the dynamics. Larger 
values of (M 1 ) and (f>( 2 K enhance the impact of cross terms 
and so magnify the deviations from the purely diffusive 
scenario. 

While cM 1 ) and <p^ indirectly quantify the amount of 
material that was introduced into the Petri dish, they 
are also probably sensitive to the microscopic shape and 
steric properties of individual ink molecules, suppos- 
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cdly distinct for cyan and magenta samples. Moreover, 
</>W and (/>( 2 ) may encapsulate three-dimensional effects, 
which although assumed to be negligible in our treat- 
ment, are responsible for a slight deformation of the drops 
along the vertical, z direction, an effect which, a priori, 
is different for the two drops. Based on these consid- 
erations, and despite the fact that care has been taken 
to keep the two drops as similar as possible at the time 
of injection into the dish, one should allow for different 
values of the unknown constants t^ 1 ) and c/)( 2 ). 

The result of the fitting procedure is reported in Fig. m 
The fitting is first performed by only adjusting the av- 
erage number densities ft- 1 ' and 4>^; here cyan is the 
molecular species 1 and magenta the molecular species 
2. The diffusivitics, D c and D m , as well as the thresh- 
olds, 7 C and 7 m , are set to the values determined above 
when analyzing the single drop dynamics (see captions of 
Figs. [3] and The time series for the widths of the two 
drops (a total of four independent curves) are hence si- 
multaneously fitted by using only two control parameters. 
The best fit values are = 0.3406 and <^ (2) = 0.2126. 

The dashed lines in Fig. [7] are the result of a generalized 
fitting strategy, where the thresholds 7 m and -f c are also 
allowed to vary. In this latter case, the fitted parameters 
result in c£« = 0.32815, (2) = 0.2233, 7m = 0.0727 
and 7 C = 0.0608. In the case of the cyan drop, shown 
in the lower panel of Fig. [7J the dashed lines obtained 
from the generalized fitting strategy which uses four free 
parameters are not distinguishable from the solid lines, 
obtained when only two free parameters are used. 

In summary, by adjusting two fitting parameters we are 
able to approximate four curves, representing the time 
evolution of the widths of the two drops along the refer- 
ence directions X and Y, with good accuracy. The slight 
deviation observed for the magenta drop is partially re- 
moved by allowing for the thresholds to change by a few 
percent, with respect to the values determined for the 
single drop case. The fact that 7 m is modified by the 




FIG. 5: (Color online) The two drops in the Petri dish, some 
time after injection. 



generalized fitting strategy has its origins in a cross-talk 
between color channels. The cyan drop induces a small 
offset in the magenta channel, when the two drops are 
present simultaneously in the analyzed image. This effect 
is responsible for the small change in j m , as compared to 
the value obtained for an isolated magenta drop. 

The results obtained point to the validity of the model 
of multispecies diffusion described in Sec. [TTJ The model 
is capable of mimicking the mutual interference between 
different chemical compounds confined to the same spa- 
tial reservoir, under crowded conditions, as seen in the 
experiments reported here. 



VI. CONCLUSION 

The problem of multispecies diffusion is of great im- 
portance in many areas of science. Particularly crucial is 
the field of biology, where distinct molecular constituents, 
organized in families, populate (for instance) the cellular 
medium under crowded conditions. Among the models 
that have been proposed in the literature to analyze diffu- 
sion in a multi-component system, the one derived in (l3j 
moves from a microscopic, hence inherently stochastic 
formulation of the problem, to macroscopic equations 
which incorporate the consequences of crowding. It does 
this by accounting for the effect of limited spatial re- 
sources, that seed an indirect competition between the 
different molecular species. The model of [TH therefore 
comprises of a set of (deterministic) differential equations 
which arc fully justified from first principles, and which 
are also expected to apply in situations when large densi- 
ties are present. In this paper we have tested this model 
experimentally by investigating the diffusion of two ink 
drops which are simultaneously evolving in a container 





<I)J5„ (L)40» (t)«l G)«l = 




(k) 70k (m)B5s (d)II(Is (o)120a 




(pj IS0J <q]145s (r) 196* (sj 245* (t) 295 s 



FIG. 6: (Color online) A collection of images displaying the 
dynamics of the two drops at different times. The images 
cover a window of observation of 295s, with an acquisition 
step of 5s. A movie of the time evolution of the two drops is 
annexed as supplementary material. 
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FIG. 7: (Color online) Upper panel: The diameters of the 
magenta drop are displayed versus time. Symbols refer to 
the experiments: open circles for the diameter of the drop in 
the Y direction and squares for the diameters along X. The 
solid lines represent the numerical fit based on the generalized 
diffusion Eqs. ©. The dashed lines are the result of a gen- 
eralized fitting strategy, described in the main text. Lower 
panel: As for the upper panel, but for the cyan drop. The 
simulations are performed over a square of linear size L, with 
L = 150. 



ter an initial evolution, which can be explained using the 
standard theory of diffusion, the fronts of the drops col- 
lide and a barrier is established at the interface which 
prevents mixing of the ink drops. This observation can- 
not be described by normal diffusion theory. The symme- 
try between the X and Y directions (defined respectively 
to be the direction of the impact and that orthogonal to 
it) is in fact broken, an intriguing dynamical feature that 
we interpret as the macroscopic consequence of the com- 
petition for the microscopic spatial resources that takes 
place at the frontier between the two drops. To test the 
model proposed in [l3j we have measured the widths of 
the drops, along the two independent X and Y directions 
and compared the experimental data points obtained to 
the numerical solution of the governing differential equa- 
tions ([5]) . The numerical time series are compared to the 
experiments by performing a constrained fitting, which 
employs a limited number of free parameters. The fitting 
algorithm incorporates information from an independent 
and preliminary analysis of the evolution of isolated ink 
drops in a similar environment. The agreement between 
theory and experiments is good, pointing to the validity 
of the mathematical model derived in [13| . We therefore 
conclude that this model accurately captures the dynam- 
ics of two, simultaneously diffusing, granular species, in 
the non-trivial situation where crowding occurs. 

The results of the experiments reported in this paper, 
together with the derivation of the model and its numeri- 
cal and analytical study, lead us to propose that it should 
be adopted in reaction-diffusion problems, where species 
are known to be spatially dense, so abandoning the in- 
appropriate Fickcan approximation. This allows novel 
insights to be gained into several important mechanisms, 
for example Turing instabilities (27j and the paradigmatic 
approach to pattern formation in biology [28[ . In [29| 
it is shown that, due to non-Fickean cross-diffusion terms, 
the Turing instability also sets in when the activator dif- 
fuses faster then the inhibitor, at odds with the con- 
ventional, currently accepted, picture. We expect many 
other applications of this microscopically motivated ap- 
proach to crowding will be explored in the future. 



filled with water. 

The experimental set up was calibrated with reference 
to the single drop case study, a simplified scenario that 
enabled us to measure the characteristic diffusion coeffi- 
cients of the ink. When instead two drops arc introduced 
into the container, a curious phenomenon takes place: af- 
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